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Abstract 

Fox's H-function provide a unified and elegant framework to tackle several physical 
phenomena. We solve the space fractional diffusion equation on the real line equipped with 
a delta distribution initial condition and identify the corresponding H-function by studying 
the small x expansion of the solution. The asymptotic expansions near zero and infinity 
are expressed, for rational values of the index a, in terms of a finite series of generalized 
hypergeometric functions. In x-space, the a = 1 stable law is also derived by solving 
the anomalous diffusion equation with an appropriately chosen infinitesimal generator 
for time translations. We propose a new classification scheme of stable laws according 
to which a stable law is now characterized by a generating probability density function. 
Knowing this elementary probability density function and bearing in mind the infinitely 
divisible property we can reconstruct the corresponding stable law. Finally, using the 
asymptotic behavior of H-function in terms of hypergeometric functions we can compute 
closed expressions for the probability density functions depending on their parameters 
a, (3, c, t. Known cases are then reproduced and new probability density functions are 
presented. 
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1 Introduction 

The notion of stable distributions was first introduced by Levy [1] in the study of Generalized 
Central Limit Theorem, and there is a nice early account of the theory in [2]. A stable law is a 
direct generalization of the Gaussian distribution and in fact includes the Gaussian as a limiting 
case. The main difference between the stable and the Gaussian distributions is that the tails of 
the stable density are heavier than those of the Gaussian density. This characteristic is one of 
the main reasons why stable laws are suitable for modelling a plethora of phenomena such as 
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laser cooling [3], turbulence [4], dynamical systems [5], statistical mechanics, signal processing 
[6], biology [7] and mathematical finance [8]. 

The goal of the present paper is to give analytic expressions for the probability density 
functions (p.d.f.'s) of the stable laws, in terms of known functions. These were lacking from 
the literature apart from a handful of well known ones. To accomplish our task we organize 
the paper as follows: 

In Section 2 we briefly review the definition of infinite divisible laws in terms of characteristic 
functions (or the Fourier transform of the probability measure) and limit our investigation to 
the subclass of stable laws. We give the characteristic exponent for the general one-dimensional 
case and comment on the role and the essential properties of the parameters which are involved. 

In Section 3 we exploit the definition of the Fox's //-function as a Mellin-Barnes path 
integral and performing the integration on the appropriate contour. This allows us to write 
its asymptotic expansion at zero and infinity (see expressions (16) and (19)). The reason for 
choosing the specific values m — n — 1 and p = q = 2 for the iZ-function is its close relation 
to the solution of the free space fractional diffusion equation. Using the Gauss's multiplication 
formula as well as standard properties of the gamma function we resum the series for ratio- 
nal values of the index a and produce general closed expressions containing the generalized 
hypergeometric functions. These expressions can be manipulated, for different values of the 
parameters, using a simple computer program running under Maple software. 

In Section 4 we solve the anomalous spatial diffusion equation on the real line with a 
fractional Laplacian consisting of Weyl derivatives and a Dirac delta distribution as initial 
condition. In this way we determine the most general form of the p.d.f. for a stable law. We 
establish the connection with the corresponding Fox function and find its asymptotics. Finally, 
we express them as finite sums of generalized hypergeometric functions. The fundamental 
solution (Green function) for the Cauchy problem of the space-time fractional diffusion equation 
as well as its relation to the Mejer G-functions have been studied by [9, 10]. 

In Section 5 we demonstrate that the diffusion equation has as infinitesimal generator of 
time translations the spatial derivative of the convolution of two terms, which under Fourier 
transformation they reproduce the characteristic function of the a = 1 stable law. The cor- 
responding integral cannot be performed exactly, unless (3 = in which case we recover the 
shifted Cauchy p.d.f. Only in the small x,/3 regime one can provide a triple series expansion. 

In Section 6 we establish a new way to classify stable laws by exploiting their infinitely 
divisible property. It is possible to write a formula that determines the p.f.d. of the stable law 
as an infinite limit of the m-fold convolution of a generating p.d.f. This expression although it 
gives new insight, from calculational view point is cumbersome due to its complexity. 

In Section 7 we present a sample of our results in the subdiffusion regime (a < 1) while a 
takes values on the Farey series 1 T n of order n — 5. In the super diffusion regime (a > 1) we 
recover all previously known results and also give new ones for general rational a. 

lr rhe Farey series T n of order n is the ascending series of irreducible fractions between and 1 whose 
denominators do not exceed n. Thus a = | belongs in T n if 

< p < q < n, (p,q) = l 
where (, ) denotes the highest common divisor of two integers. 
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2 Preliminaries on a-stable laws 



Consider a probability measure fi on R n and its characteristic function [11, 12, 13] 

£(p) = F\j*](p) = I e l< ^>fi(dx), p G R". (1) 
Jm. n 

A probability measure fi on R n is called infinitely divisible if 

Vm eiV, 3 // m , J"[/i m ] : ^[At](p) = (-^[Atm](p)) m 

where fi = fi m * ■ ■ ■ * fi m is the m-fold convolution of fi m with itself. If the measure fi is infinitely 
divisible then there exists a unique continuous function i/j : R™ — > C, called the characteristic 
exponent of fi, such that ^(0) = and 

F[fi]{p) = e* (p) , peR" 

The Levy-Khintchine representation or Levy-Khintchine formula states that a probability mea- 
sure on R n is infinitely divisible iff we can write the characteristic exponent in the form 

^(p) =i <p,r> -- <p,Ap> + ( e i<p ' x> - 1 -% < p,x > l| x |<i) u{dx) (2) 

where r G R n , A is a symmetric nonnegative-definite nx n matrix, called the Gaussian covari- 
ance matrix, and v is a a-finite Borel measure on Rq := R n /{0}, called the Levy measure, such 
that 

I m in { l,|W|> W <oo. (3) 

The triplet [r, A, z/] is unique and will be called the generating triplet of the infinitely divisible 
probability measure fi. If A = then fi is said to be purely non-Gaussian. 

A subclass of infinitely divisible laws is the stable laws class. Suppose that X, Xi, ■ ■ • , X m de- 
note mutually independent random variables with a common distribution F and S m = Y^hLx 
The distribution F is stable if for each m G N there exist constants c m > and r m G R such 
that 

S m = c m X + r m (4) 

and F is not concentrated at one point. F is stable in the strict sense if r m = 0. The symbol 
= means that the distributions of S m and X are identical up to scale and location parameters. 
The norming constants are of the form c m = with < a < 2 and the constant a is called 
characteristic exponent of F or index of the stable law. 

Let < a < 2 and fi be an infinitely divisible and non-trivial on R n probability measure 
with generating triplet [r, A, v\. If fi is a-stable then there is a finite non-zero measure A on 
the unit sphere S = {x £ I" : |x| = 1} such that 

{%) A = and v{B) = J s \(d£) / °° l s (r£)^ for B G B(R n ). 

(ii) fi(p) = exp [— J s \ <Vii > ~~ ^tan(^))sgn < p, £ >)A(d£) + i < r, p >] for a ^ 1 and 
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(in) fi(p) = exp [-J s (\ < P,£ > | +ijf < P,£ > In | < p,£ > |)A(d£) + i < r,p >] for 
and r Gl n 



ft 



In the one-dimensional case (n = 1) one can prove that the characteristic exponent has the 
form 

,fs ■ f c\p\ a [l-i(3 sgn(p)tan(f )] if a ^1,2 _ 

^(p)=«rp-^ ' ' v\ „ 2 J/j, n /uni ;* ' i ( 5 ) 



c|p| [l + i/3f sgn(p)ln(|p|)] if 



ft 



where 



P = C+ , ° » and C= 77F7T^ — r • La\ ( c + + c ~)' for a ^ 2 ( 6 ) 
c+ + c_ 2I\1 + a) sm (^) 

with c + , c_ > and c + + c_ > 2 . Note that when a = 2 then = 0. 

The collection of the four parameters (a, (3, c, r) is called the stable law parameters and 
completely determines the distribution as follows: 

Characteristic exponent a. This parameter determines the degree of leptokurtosis and the 
fatness of the tails. For a stable real-valued random variable X it can be shown that 

E\X\ < oo 1 < a < 2. 

When cc < 1 the means becomes infinite. The variance for a G (0, 2) becomes infinite or unde- 
fined while all moments of a random variable X become finite iff a = 2. Also the moments of 
order less than a are positive and have a finite limit, namely E\X\ k < oo, < k < a 3 . 
Skewness parameter (3. This parameter characterizes the degree of asymmetry of the Levy 
measure and takes values in the interval [—1, 1]. The measure v is called symmetric if (3 = 
(or c + = c_) and the a-stable distribution is called stable a-symmetric 4 . 
S'ca/e parameter c. This parameter ranges into the interval (0, oo) and measures scale in place 
of standard deviation. 

Location parameter r . This parameter saturates the set of real numbers and shifts the distri- 
bution to the left or right. If 1 < a < 2 then r equals to the mean of \x. When < a < 1, 
although, the mean is infinite it serves as an index of the location of the peak of the stable 
distribution and is identical to the drift of \i. 

A stable law generated by (a, {3, c, r) is often denoted by S a (P,c,r). In the present work 
our law will be first generated by S a (c) and then we will study the most general case. 

The parameter space (a, b) of stable p.d.f.'s with the centering constant b restricted in the 
region [15] 

1 7 1 r cx , < a < 1 /_s 

l^{ 2 - a , l<a<2 < 7 > 

is depicted in the following figure 

2 The same result can be recovered if one uses the absolutely continuous Levy measure 

u(dx) = (c + l x>0 + c-l x<0 ) \x\~ 1 ~ a dx. 

3 It is also true that a symmetric a-stable random variable has finite negative-order moments — 1 < k < 
[14]. 

4 In general a measure fj, is symmetric when fi(B) = fj,(—B) for B € S(R"). In n = 1 the rotation invariance 
is tantamount to symmetry. 
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Figure 1: The parameter space of all stable p.d.f's on the r = plane. The points of the axis 
/3 = represent the stable a-symmetric p.d.f.'s with the property / Qj/ g=o(^) = f a ,/3=o(— On 
this axis are located the familiar Cauchy (a = 1), Holtsmark (a = |) and Gaussian (a = 2) 
distributions. 



3 Fox's //-function and generalized hypergeometric func- 
tions for the law S a {c = K a ) with rational a 

Fox [16, 17] defined the if- function in his studies of symmetrical Fourier kernels as the Mellin- 
Barnes path integral 



(a 1 ,Ai),(a 2 ,A 2 ),- ■■ ,(a P ,A p ) 
(b 1 ,B 1 ),(b 2 ,B 2 ),---,(b q ,B q ) 



1 

2m 



X{s)z s ds 



(8) 



where the integral density x( s ) is given by 

Yl™ 1 r(b i -B i s)nU r ( 1 - a i + A i s ) 



X(s) 



U!= m+ i r(i -h + b iS ) nu+i T ^ ~ A ^ 

m, n, p, q are integers satisfying 

< n < p, < m < q, 
Bi, Aj are positive numbers and bi, cij are complex numbers such that 



(9) 



A j {b h + v)^B h {a j -l-X), u,X = 0,l- 



h 



m; j 



n. 



This condition implies that the poles of F(bi — B^s) and T(l — a,j + Ajs) form two disjoint 
sets. C is a contour in the complex s-plane which runs from s = oo — ik to s = oo + ik with 



k > 



| Im6j| 



j 



n and which encloses the poles 
(k + ") ■ , 



m 



but none of the poles 



(dj — 1 — v) 

a; 



3 = 1, 



makes sense and defines an analytic function of z in the following two cases: 

(0 if 

q p 

M = B i ~ Yl A 'J > °' V2 ^ 
i=l j=l 



(10) 



(zz) If 



M = and < \z\ < R with # = J^(5i) Bi jJ(A,•)- 



A, i 



i=l 



Schneider [18] has introduced into physics the Fox's iJ-function as analytic representations 
for the Levy distributions in x-space, and as solutions of the fractional diffusion equation 5 . If 
one tries to solve the one-dimensional anomalous diffusion equation equipped with the following 
initial and boundary conditions 6 



df(x,t) 
dt 



K a A{a)f(x, t), iei,i>0,ae(0,2) 
5(x) 



, lim f(x,t) = 

\x\— >oo 



(12) 



then the solution (or propagator) reads 

1 



f{x,t;a) 



a\x\ 



H. 



1,1 

2.2 



X 



1 i ) 

(1,1), 



(13) 



which expresses the a-symmetric stable p.d.f. in terms of the Fox's i?-function. Thus our study 
will be focused on the i?-function 



H(z) 



H. 



i.i 

2.2 



:i,§) 



r(i- a )r( 



2m j c r(|)r(i 



j-z s ds, z 



\x\ 



(K a ty. 



(14) 



The simple poles of T(^) and T(l — s) are given by the disjoint sets of points 

P(s) = {s u = —au, v — 0, 1, ■ ■ •} 

Q(s) = {s„ = 1 + v, i/ = 0, 1, - - ■} 
We distinguish the following two cases: 



5 In nature there is a diversity of diffusion processes for which the fc-moment of the displacement grows not 
linearly with time but follows a power-law pattern of the form E\X(t)\ k ~ , < k < a. 

6 In equation (12) K a is the generalized diffusion constant having dimensions [L] a [T] -1 and A is the operator 
(33) with (3 = and r = 0. This is the generalization of the Laplacian to a fractional order. 



6 



(i) Asymptotic expansion of H(z) near the point z = oo. Applying the residue theorem 
clockwise we find 



H(z) = f2Res{ X (s)z s ;s m e P(s)} + ^- [ X (s)z s ds-^-[ X (s)z s ds 

2m Jci 2m Jc 2 



(15) 



where the contour C has been replaced by the rectilinear contours C\,Ci running from 
a ^ a + ik ^ — oo + ik and o ^ o — ik ^ — oo — ik respectively. It can be proved that 
the error terms on the right-hand side of (15) vanish. In this case we obtain the algebraic 
asymptotic expansion 



H. 



1.1 

2,2 



(hi) 



^' i' 



a 

7i 



m=l 



_ x r(l + ma) . ,ir 
— ; — sin —ma)z 



(16) 



It is worth noting that all the points of P(s) contribute in the large z limit. For the series 
(16), applying the ratio test, one has 



Pi 



lim 

m— >oo 



r((m+ l)a) sin(f (m + l)a) 



< lim 


T((m + l)a) 


mT{ma) 


m—*oo 



if < a < 1 

1 if a = 1 (17) 
oo if 1 < a < 2. 



mT{ma) 

Thus the series converges absolutely for every value of z ^ in the interval 

(-Ri, Ri) = (-oo, oo) if < a < 1. (18) 

(u) Asymptotic expansion of H(z) near the point z = 0. The function if(z) is analytic in the 
s-plane for a G (1, 2] since then M = 1 - ^ > 0, 7^ 0. Also for a — 1, M — 0, and if 
is analytic for < \z\ < 1. In this case we find 



H. 



1,1 
2,2 



(1,1) 
(1,1) (1,1) 



'2m+l' 



Y: Res{ x (s)z- s m g Q( a )} = - E(- 1 ) m f7j^TT 



2 



2m+l 



777=1 

oo 

-£< 

7T 

771=1 



771=0 



T~Yl -I- — ^ 
■l)™- 1 a sin(-m)^ 

7 ri + m V 2 ; 



(19) 



In this case the even numbers of Q(s) give a vanishing result and 

r 



P2 



lim 

777— > OO 



■ 2777+3 ' 

a ' 



00 if < a < 1 

1 if a = 1 
if 1 < a < 2. 



T(2m + 3) 

The series converges absolutely for every value of z in the intervals 

(-1,1) if a=l 
-R 2: R 2 )-i ( _ QO)0o) if 1<a < 2 . 



(20) 



(21) 



We assume now that a is a positive rational number thus it can be written as: a = 2 , p, g G Z + 
(the symbols p, g should not be confused with those used in the definition of the if-function). 
Relation (16), using the substitution m — nq + I, I — 0, • • • , q — 1, can be casted into the form 



H(z) 




a x ^ 

7T ^-^ 



»r(np+f + l) * 
^ -^7 —7—7^ sm( — ( n 



77=0 



T(nq + l + V 



)p)e l7Tnq z- np 



(22) 



Using the multiplication theorem of Gauss one can write the ratio of the gamma functions as 



r^ + / + i) YAW re:Jr(^fi) nS^)n 

where (a) n is the Pochhammer's symbol (a) n = r ^"- > and 

fc + 1 Z , Z+s+1 

a fc = + -, 6 S = . 24 

p q q 

Combining (22) with (23) and assuming that p is odd, we obtain 

H(z) = al (2tt) V 1 X: n f° r(afc) ^in(^) f] ^e^S)' 



2m 



00 TTP~l/'„ \ /„,P„mt(<J+S)\ 2m + 1 - 



+ e-^cos^sin^VlS^^f^^) ) 



IE(M S 

27r 2 \ — 2L-^ sin +e 2 ^cos sin — p 

1 ; t^v^; rc=;r(6 s )V l 2g ; 



9-1, 



The same result holds for p even and the only difference is the swapping of the two terms in 
the parentheses. Of course the hypergeometric series converges absolutely when < a < 1 — ^ 
which implies that a 6 (0, 1). 

For the case 1 < a < 2, following similar steps, we can write the ratio of gamma functions 

as 



(26) 



r (^) (pV (IX ( (2g)^ \ n£ 1 r(a.) Ull^n 

r(2m + i) a^ 1 ^ 1 -^ W W V(2p) 2p V nST'rfe) nX 1 ^). 

where = ^(2/ + 1) + ^ and b s = ^(2/ + 1) + ^. The relation (19) then becomes 

{ ' y ] p > \p p2p - q ) nS rfe) 

x 2 q F 2p -i (l, a 1; ■ • • , a 29 _ i; 6 1} ■ • • , 6 2p _ i; M^ e -^ . (27) 

The hypergeometric series 2< ? -Z 7 2 P -i converges absolutely if 2q < 2p — 1 or equivalently when 

a > 1 + — , geZ + (28) 

which implies that a G (1,2]. It also converges when 2q = {2p — l) + l=^a = l provided that 
\z\ < 1. 



7 The multiplication theorem of Gauss states that r(mz)(2 7 r) (m - 1 )/ 2 = m mx -^ 2 T(z)T(z + ±)---T(z + 
s=i), Vm e TV. 



4 Fox's H-function for the law S a ([3,c = K a ,T), ol G 
(0, 2] and a/1 

We consider the generalized one-dimensional anomalous and anisotropic diffusion problem with 
the following initial and boundary conditions 



df(x,t) 



dt cos 



(K-^V: + K+ x V^f(x,t)-r dJ 



dx 

x £ R, t £ (0, oo), \imf(x,t)=5(x) 



A{a,P,r)f(x,t), 
lim f(x,t) =0(29) 

\x\— >oo 



where -RT are diffusion constants satisfying AT > 0, K + + K > and the dispersion term is 
proportional to the constant r having dimensions [L]/[T] 8 . Also by definition [23, 24] 



clef 



T(m — a) \dx 







f(y,t) 



(x-y) 



a—m+ 



(y-x) a ~ m+1 



are the nonlocal fractional left-handed (right-handed) Weyl derivatives. In (30) m = [a] + 1 
with [a] representing the integral part of a. Note that when a is an even integer then the two 
derivatives are localized and equal while for odd integer values of a both derivatives appear 
opposite in signs. 

The Fourier transform of A(a, j3)f(x, t) for fixed t £ (0, oo) and / £ 5(R) 9 is given by (see 
proof at Appendix A) 



A{a,f3)f {q,t)= -A' Q |g| a (l-z/5sign(g)tan(^— 



%rq 



f(q,t) 



where 



K + — K~ 

K a = K + + K~>0, p= : — £ [-1,1], r£ 



K~ + K+ 

The operator A(a,f3) could be casted into the equivalent form 

K„ 



(31) 



(32) 



A(a, 0) = ((1 - ;) ,K + + X V«J - r±. 



(33) 



The f(q,t) satisfies the Fourier transformed initial value problem 
df(q,t) 



dt 
/(?,0) 



A a |g|° ( 1 - i/5sign(g) tan j 



with solution 



f(q,t) = e - K ^\ at { 1 - i ^ n ^ tan ( 2 T)) +iTqt . 



(34) 



(35) 



8 One might wonder if it is legitimate to add a Laplacian term to the operator A(a,f3). This suggestion is 
prohibited by the fact that we consider only a-stable laws. 

9 <S(R) is the set of all infinitely differentiable and rapidly decreasing functions on R, namely 
sup xG r\x" {D m f)(x)\ < oo, Vto, n = 0, 1, ■ • -. This space is usually called the Schwartz space. 
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The propagator is thus given by 

f(x,t) = ^[ftM = — r e^-^e-^F^^M^&Vdq. (36) 

27r J-oo 

We will first study the small x — rt expansion. We expand cos(ga;), sin(gx) in finite Taylor series 

(— q 2 x 2 ) n g(qx)(qx) 2m+l 



cos 



(2n)! (2m + 1] 



W " ^ (2n + l)! + (2m + 2)! (3?) 

where, by the generalized mean value theorem, the functions g(qx), h(qx) are bounded by the ex- 
treme values of the (2m + l)th ((2m + 2)th) derivative of the cosine (sine), thus \g(qx)\, \h(qx)\ < 
1. By substituting these series into (36), integrating term-by-term and taking the limit m — > oo 
we obtain the complete asymptotic expansion 10 



>"(l + 7 2 )^^ r(n + l) V 



f{x, t- a, /5, r) = — — > J-IT' 1 '^: ', «' sin ) z" (38) 

7TZ (-ftT a t) ' 



where 

/ OJ7T \ 2 X — 7~£ 

7 = /3tan( — ), 5 = 1-1 arctan7, and z = j j-. (39) 

V 2 J air {K a t)« (1 + 7 2 )^ 

It is evident from (39) that 5 G [0, 2] and can acquire the three integer values {0, 1, 2} provided 
that j3 G { — 1,0, 1}. The p.d.f. (38) can be reproduced by the following H-function 



H(z) = H$ 



(1,1) (i,|; 



rn-s)r(i) sj \x-rt\ 

' Ka z s ds, z = l - 1 (40) 



2™ Jc r(f )r(i - f ) (jr a *)S(i + r») 

with asymptotic expansion near the point z = given by 

ffW = ^(-i rl -X±|isin(^),", l< a <2. (41) 

n=l ^ ' 

In contrast to (19) all the points of Q(s) now contribute. The asymptotic expansion of H(z) 
near the point z = oo is 

^)^EHr|^-(™*)W- <>«,<!. (42) 

n=0 v 7 



10 Thc derivation of this expression is based on the integral formula 

f x^ 1 e~ ux suv(ax)dx = — — sin f/i arctan — ) , Re/i > —1, Reiv > |Ima| 

Jo (v 2 +a 2 )2 V vJ 

and a similar result for the cosine, provided that Rc/i > 0, Rc^ > |Ima|. 
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For the symmetric /3 — case we recover expressions (16) and (19). When a is rational we can 
also express (41) and (42) in terms of hypergeometric functions as follows 



H(z) = -JV)-- 1 ^ 



p—l/3. 



qte i* z \ y\1 = \ T(a k 



1=0 



V 



u p s =lm 



l 2 ( 18 1 ) g F p _i I 1, ai, ■ ■ ■ , a q -i, fei,---, bp-i, 



e i ^ {l5+1 \F p _ 1 ( l,ai,---,a,_i;6i,---,6p_i; 



; ivrp(l-|)^q 2 p 



1 < a < 2 (43) 



with a fc = ~ + ~ and o s = ^ + ~. Similarly for the large z expansion we obtain 



H(z) 



=0 



P+i^g | 1, «o, • • • , a P -i; oo, • • • , bg-i, 

e t7r(g-^y 

g<?|z|P 



e -if(^«+l) 



3+1-^9 ( 1) a 0) ' ' ' 5 a p-l] bo, - • ■ , 9 _i 



9-1) 



< a < 1 (44) 



where a fc , o s are given by (24). 

5 The law S a=1 (/3, c = K, r) 

The operator we consider in this case is 



A(a = 1, (3)f(x, t) = ~ {K x {s + (3h) * f{t)){x)) 



where 



s(x) = —z-, h(x) 



1 t 1 



+ 2C 1 8{x) ) j . 



2ttV v ' 2vr 2 
The first convolution is the Hilbert transform of the function / defined by 

fOO 



Hf(x,t) = (s*f(t))(x) 



2tt 2 



f(y,t)s(x -y)dt 



(45) 



(46) 



(47) 



with P.V. representing the Cauchy principal value of the integral. The Fourier transform 11 of 
(47) is 



T[Hf]{x) = F[s](x)T[f]{x) = —sign(q)f(q,t). 



(48) 



L We have absorbed the coefficients of the Fourier transform into the exponential thus defining 



F[f\ = f(p) = / e 2i ^ x f(x)dx, T-^f] = f(x) 



f(p)dp. 
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The second convolution term in (45) has Fourier transform (see Appendix B for the proof) 

F[h*f\(x)=F[h](x)r\f\(*) = -^H\q\)f(q,t). (49) 
Hence, the initial value problem is equivalent to the Fourier transformed 
df{q, t) 



—K\q\ ^1 + i/3— sign(g) ln(|g|)^ + liixqr 



dt 

f(q,0) = 1 (50) 

with solution 

f(q,t) = e - X l 9 l*( 1+ ^f si g n (9) ln (kl))+2*vrgrt_ ^ 

The propagator is then given by the absolutely convergent integral 

POO 

f(x,t)= / e - 2i ^ x - T t) e -K\ q \t(i+i^sign iq )HW) dq _ (52 ) 



OG 



If we set (3 = in (52) we recover the shifted Cauchy p.f.d. 

2K x t 

f (X ' t] ~ {K.tf + A^{x - rtf (53) 

6 An alternative way of classifying stable laws 

The symbol of the operator (33) is given by 

77(g) = -K a \q\ a (1 - isign(g)7) + irq. (54) 

From (54) the real part Re(?7(g)) < 0, Vg E E, thus we define h x : R — > C, A > by 

f°° 1 

h x (q) = £[e tri{q) ](\) = / e- tx e^dt = - (55) 

Jo A - 77(g) 

which is positive definite. The mapping q — ■> h\(q) is continuous and applying Bochner's 
theorem there exists a finite measure on £>(M) such that 

hx(q) = /£(?) = ^ r e ~ iqX ^)- ( 56 ) 

Hence, 

POO 

fi x (x) = / e^h x {q)dq. (57) 



It can be shown that the operator A(a, (3, r) is the infinitesimal generator of a strongly contin- 
uous semigroup of operators and its resolvent 1Z\(A) satisfies [13] 

K x (A)ip = pi X *ip, (58) 

12 



or equivalently 



1 ' 1 



I - -A J i) = A/i A * 4>. (59) 



By applying on the righthand side of (59) m-times the operator (^1 — jA^j , setting A 
and taking the limit m^oowe have 

lim 1 1 -—A) V = e~ tA ip = lim I (^] n m/t * • ■ ■ * /i m/t ] * ^ 

m^oo y m J m^oo I V I / !✓ / 



f 



m— times 

= /*V (60) 

where / is the p.d.f. of the stable law S a (/3, K a ,r). Thus, a stable law can be determined by 
fix instead of /, using relations (55) and (57). As an example we consider the law S^O, K 2l 0). 
A simple calculation gives the p.d.f. 

u A (x) = — ^=e~^ x , x>0. (61) 

Although (57) seems to be elegant it is unattractive for calculations since the corresponding 
convergent integral only exceptionally gives a closed expression. 



7 Known and unknown results 

Expressions (25), (27), (43) and (44) can be used to calculate p.d.f.'s for arbitrary rational 
values of a G (0, 2] in terms of generalized hypergeometric functions. In particular, a sample of 
stable symmetric p.d.f.'s is given in Table 1 for the subdiffusion regime with a G (0, 1) and in 
Table 2 for the super diffusion regime with a G (1,2]. Previous known results are reproduced 
and new ones are also presented. 



(a) The law 5 1/2 (0, K 1/2 , 0) 

Let us first consider the particular case (p, q) 
From equation (25) one can easily show that 



^1, 2) that gives the index value a = 1/2. 



H(z) 



1 1 / e i7r \ 1 ( °° 

1 



471Z 

8-Jttz 



COS( 4^ )lF2( 4 ; 2'4' 64.- 
4^ oFl( 2 ; "64^ ) + 0jPl( 2 ; 



sin(f(n+f)) 1 

r( n + ^) W) n 

1 , 1 



sin l — 

12z Az' 



QAz 2 ' 



1 



2y/2ir\z 
1 



cos(^)C(— = 
Az ! K ^hTz 



sin( — )S( . 

Az ! K ^hTz 



4V2 



irz 



cos(l) + sin(l) 



F( 3 - 3 7 



4' 2' 4' Q4z 2 



(62) 
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Where C(x), S(x) are the cosine and sine Fresnel integrals given by 



C(x) 



S{x) 



J ™s(-t)dt = x 1 F 2 (-;-,-,--x) 

/ sm(-t 2 )dt = -x 3 1 F 2 (-:-,-, x 4 ) 

J v 2 ' 6 iV 4'2'4' 16 ; 



x tM; /V^ 1 "*) -4 oFi(l,-^xH)dt (63) 



r(f ) 2 7 v 1 n 4' 16' 



which are odd functions of x. The function f(z) = -^H(z) is indeed a p.d.f. since it is 
positive definite Vz G (— oo, oo), integrable and satisfies 

/(*)<fe = 1. (64) 

> 

The more general one-sided law 61/2(1, ^1/2, T ) is found to correspond to the p.d.f. 

/(^) = — L— e -£, z>0. (65) 

The case f3 = —1 gives a vanishing p.d.f. since 5 = 0. This result holds independently of 
the value of a. 

(/3) The Cauchy law S^O, K x , 0) 

Setting in (27), p = q, yields 



9-1 



ire*) = VoaiOiivn^-i) 



Z 

IX 

1=0 



K 1 TTl? ! ) li: '° (1;0;( - z2)S) - (66) 



When g = 1 this corresponds to the Cauchy p.d.f. 

= Jl£T(z) = !— |— , zG (-00,00). (67) 

az 7T (1 + 2 J 

with applications , e.g. in spectroscopy. For the general law Si(Q, Ki,t) the p.d.f. is 
given by (53). 

(7) The Holtsmark law 5 3/2 (0, K 3/2 , 0) 

In this case we set (p, q) = (3, 2) in (27) and arrive at 

*2§ * ^,/2 3vnLir(M2/+i) + !) 



*w - ^ EC-*) 1 ^ 



v/3fe V27 ; nLir(|(2/+i) 



1 11 31 11 5 

L,-(2/ + l) + i ,---,-(2/ + l) + i ;-(2/ + l) + -,---,-(2/ + l) + - 
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6 6 

z 

TV 



e 3 ™z 6 



r ( 2) 2 F,(A,H ; I,I,5 ; _A z6) 

V v 12 12 3 2 6 729 ; 



3^4(1, 



2 3 ^ 4 ^' 4' 4' 6' 3' 3' 6' 729 



z 6 ) 



14^3 7T , n ,19 13 3 7 5 



486 r(|) 



12' 12' 2' 6 3 729 







The corresponding p.d.f. was discovered in physics by the Danish astronomer Holtsmark 
in 1919 [26, 27]. It was the outcome of his efforts to study the stationary distribution 
of the force acting on a star, per unit mass, due to the gravitational attraction of the 
neighboring stars. 

(5) The Gaussian law £2(1, K 2 , r) 

Substituting p = 2q in (27) we get 



H(z) 



z(27l 



1 9-1 



2q-l 



9N I T-r2g-l r / 2(H-fc)+l 



z=o 



n ^-i r (^±i) 



X 



(4g)4g 



z(2tt 



,3-1 ^ 



V 7 ? 



Em)' 



z=o 



v2\ ' 



2g-i p/Z+s+i/ 



*J IIS r( 



2<z 



2<H 



-^4(2-1 I 1, Ol, " " " , CL2q-l', bl, • • • , &2p-l! 



r 4, 



(4 3 g 



2W 



where a fc = ^1 + A anc i fr s = 2L±1 _|_ jl_ Setting g = 1 in (69) we have 



2p 1 2p • 



(69) 



tf(z) = — 

'71 



oxtail) ~ ~T° Fl ^2 ] M J 



e * . 

7T 



—j= ( cosh(— ) — sinh(— ] 



(70) 



The associated p.d.f. is 



and changing variable to z 



/(*) 



2A 



e 4 



^t= we recover the traditional one 



/OM) 



y/AnK 2 t 
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Table 1: The Fox's if-function for the Farey series T n of order n = 5 excluding the first ° and the 
last \ member of the series. The functions csc(z) and I v {z) are the cosecant of z and the modified 



Bessel function of the first kind given by: I v {z) = e lw ^J u 



i 



i>+i) 



a 



H(z) = a\z\f(z) 



1 csc (f ) sin (ik) 



H^p (1 3 2 1 7 4 9 . 



"0^7^55 io> 5> 2» 10' 5' 10' 25-10 8 z 2 



25 2^r(|) 

1 csc( ^) sin(f) /_3 2 1 3 4 _9_ 11. 



) 



25 



"5~~~T 0-^7^) |) f ) f ) g ) fo"> fg! ~25-10 8 z 2 - 



1 r(§)sin(§£) 



50 



3 

.Z^ T „ 

1 r(|)sin(^) 



0-^7 ( 



2 1 3 7_ _9_ 11 6. 

5' 2' 5' 10' 10' 10' 5' 



) 



1 

25-10 8 z 2 
~25-10 8 z 2 - 

...Lr) 

1 csc(§)sin(%^ , x 4 _9_ 6 13 7 3. 1 



150 



PCi 3 1 4 11 6 13. 
O^ 7 ^' 5' 10' 5' 10' 5' 10' 



1 l 2 pji. 3 7 6 13 7. 



+ 600 ztt 1 ^^] 5> 10' 5' 10' 10' 5' 10' 5' 25-10 8 z 2 - 

7? (J. 4 _9_ 6 13 7 3. 1 
O^ 71 - 10' 5' 10' 5' 10' 5' 2' 25-10 8 z 2 - 



+ 



15000 2 « r( 4 } 
1 csc(^)sin(f^ ) _ j-, /4 g 11 13 7 3 



45000 



1 



5' 10' 10' 10' 5' 2' 5' 25-10 8 z 2 ' 



1 r(|)sin(f) 



210000 



52„7? (JL 11 6 7 3 8 17. 

^ 7 m i 1 n i rr i rr i o i r i i n i 



" or7 ' 10' 10' 5' 5' 2' 5' 10' 25-10 8 z 2 



) 



1 r(|)sin(^) 



1260000 



V 7 2sin( 



"9 

7T 



F 7 ( 



116133817 9. _ 1 ' 
10' 5' 10' 2' 5' 10' 5' 25-10 8 z 2 < 



16 



*T(4) 



EZ„J? 5 (I 3 13 7. 



1 



4' 8' 2' 4' 8' 4194304z 2 



\ 1 V?) p /3 1 5 7 9._ 1 

/ 32 ^1^0^518' 2' 8' 8' 8' 4194304z 2 - 



~ 3i 3 0-* 5 V 2 ' 8 ' 4 ' 



32 "S^ 

+ 1536 57^77; o^sv 



15 3 9 5. 

2' 8' 4' 8' 4' 



1 



4194304z 2 



3 7 5 11 3. 

4' 8' 4' 8 ' 2' ' 



" 41943042 2 



\ 1__L c 1 /i . 5 3 7 9 5 11. 

> 96 ztt 1 "^ 1 ' 8'4'8'8'4' 8 ' ~ 

\ 1 yg p / 7 9 11 3 13. 
7680 1 _0-^5^s ' 8' 8 ' 2' 8 ' 



+ 



1 r(|)sin(f) 



z? 



1 A 
4194304z 2 > 

1 

4194304z 2 - 



15360 



Z^7T 



7? (2. 5 3 
0^510, 4 , 2 , 



1 



8' 4' 2' 8 ' 4' 4194304z 2 ■ 



zT?T(|) 



/l 1 5. 



1 



11664z ; 



1 Vsrd) . 

— 0^31 



2 7. 



18 



3^2' 3' 6' 



Z37T 



4_J_ J_ T? 5 7 4. 

18 ZTr 1 ^ 4 ^ 1 ' 3' 6' 6' 3' 



1 



11664z 



162 



«r(|) 



18 Z7T 

F q (5 4 3. 1 \ \ 1 r (D TP- ( 7 3 5. 

316' 3' 2' 11664z 2 ' n 



648 4 

Z 3 7T 



01 3Vfi' 2' 3' 



1 



2' 3' 11664z 2 • 



1 TCT^WTw^i Ff7 17.1 3 21 7 4 9. 

„,i>"H 5 J &m V5^ 2 - r7 ^20' 20' 5' 10' 5' 2' 10' 5' 10' 



' z §r(^) 5 Vs ' 7 ^20' 20' 5' 10' 5' 2' 10' 5' 10' 39062500z 4 - 

"25 1 C5) Sm l — J^p2^7l20, 20' TO 

Z5 7T 

j 2 csc(|)sin(-^) tp ,iA Zl z 6 1 \ 

"tiok L,.,, 2- r7 ^20' 20' 5' 2' 5' 10' 10' 10' 5' 39062500z 4 ' 



v 20 ' 20 ' 10 ' 5 ' 2 ' 5 ' 5 ' 10 ' 10 ' 39062500z 4 
•11 21. 2 1 3 7 9 11 6. 1 



125 zSr(|) 

2 r(|)sin(|) p / 13 23.1 3 X 6 13. 



"5" 

Z 5 7T 

r 27r ^ 



+ 



125 

7 csc(^) sin(|) 



T2" 



*rr(|) 



7 ^ 20' 20' 2' 5' 10' 5' 10' 5' 10' 39062500z 4 

p /17 27. X 4 

2 - r7 v20' 20' 10' 5' 

9 11 13 7 3 8. 



) 



_9_ 6 13 7 3. 1 

9375 ia r , 3 N £J - 20 ' 20 ' 10 ' 5 ' 10 ' 5 ' 10 ' 5 ' 2 ' 39062500z 4 
z & HjJ 

1 r(t^ )sin (T : ) p /19 29.4 9. 11 13 7 3 8. 1 \ 
3125 H 2 - r7 V20' 20' 5' 10' 10' 10' 5' 2' 5' 39062500z 4 > 

Z 5 7T 

CSC (f)sin(^) 



11 cscjjjsm^i p /21 31. J. 11 6 7 3 8 17. 1 
" r 656250 T^ETTIT 2 7 ^20' 20 ' 10 ' 10 ' 5 ' 5 ' 2 ' 5 ' 10 ' 39062500z 4 - 



z^r(|) 



13 r (s) sin (f )_ Z7^^23 33.11 6 13 3 8 17 9. 



1312500 



1 



2V2ttz 



— — 



20 ' 20 ' 10 ' 5 ' 10 ' 2 ' 5 ' 10 ' 5 ' 39062500z 4 



) 



+ 



1 



4V2ttz 



cos(^) +sin(^) 
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3 II 9 r (f ) si ^(ffr) p ( 4 13 23 14. 1 3 2 1 7 4 9. 729 \ 

5 25 1 4-^71 15 ) 30 ' 30 ' 15 ' 5 ' 10 ' 5 ' 2 ' 10 ' 5 ' 10 ' 15625-10 4 2 e > 

9 cgc (f )sin(-f) p /ll _8_ 13 31 . _3_ 2 1 3 4 _9_ 11. _ 729 \ 
125 6 4 4-r7l3 , 15 , !5, 30 , 10 ; 51 2' 5' 5' 10' 10' 15625-10 4 2 6 / 
4 ^' 

I 18 r (s) sin (Tg) p (7_ 19 29 17.2 1 3 _7_ _9_ 11 6. 729 \ 

~ l ~125 I 4 7 U5' 30' 30' 15' 5' 2' 5' 10' 10' 10' 5' 15625-10 4 2 6 > 

21 c ^ c (-#) sin (|) p /17 11 16 37.1 3 J_ 4 11 6 13. _ 729 \ 
625 z ¥y(±) 4 ^ 30 ' 15 ' 15 ' 30 ' 2 ' 5 ' 10 ' 5 ' 10 ' 5 ' 10 ' 15625- lO 4 ^' 
3 1 p /2 5 1 7 4. 3 X 4 _9_ 11 6 13 7. _ 729 \ 
100 i^r" 5 8 ^3' 6' x ' 6' 3' 5' 10' 5' 10' 10' 5' 10' 5' 15625-10 4 2 H ' 
117 r (s) sin (f) 77 /23 14 19 43. 7 4 9 6 13 7 3. 729 



"15625 T5 4-^7^30) 151 15; 30 ' 10' 5' 10' 5' 10' 5' 2' 15625-10 4 2 B • 

Z o T 



I 11c h 5 ;^ 5 ; p /; 
+ T5625 3 4i ^7l-- _ 

66 csc (5) sln (io) p /13 31 41 23.4 _9_ 11 13 7 3 8. _ 729 \ 
"T" 78125 z Qt(±) 4 ^ 15 ' 30' 30' 15' 5' 10' 10' 10' 5' 2' 5' 15625- 10 4 2 b / 

171 r (|) sin (T-) p /29 17 22 49. _9_ 11 6 7 3 8 17. _ 729 \ 
156250 M 4 - r7 V30' 15 ' 15 ' 30 ' 10 ' 10 ' 5 ' 5 ' 2 ' 5 ' 10 ' 15625-10 4 2 6 ' 

I 561 csc(^)sin( 1 g-) p/16 37 47 26. 11 6 13 3 8 II 729 ■ 

"""31125000 ^TTT(i) 4 ^ 15 ' 30 ' 30 ' 15 ' 10 ' 5 ' 10 ' 2 ' 5 ' 10 ' 5 ' 15625-10 4 2 B - 



I (^M^)^(^) + 9^ Sinh (2^Kl(2^) + 9^CO S h(^)/ | (^ ; 

-37 (6<«rfi(^|(3^)) + SF sinh (2^)/|(2^) + 3^ cosh(^)/ f ( 3 
+ 8iV (cosh(^)/ | (^) +9z 2 sinh(^)/ | (^) + S inh(^)/4(^ 
- SEP (^M^lfe) + IS^smhC^r)/^^) +s inh( 1 ^)/ | (^ 

„,:t, — : — / Mtt ; 



2^)) 



3 9 r (4~) sill (-f) p /_7 11 19 23.1 3 1 3 7. 729 

4 16 3 4 5 v 24 ' 24 ' 24 ' 24 ' 4 ' 8 ' 2 ' 4 ' 8 ' 2621442 B 

9 V2 r/1111 13.3157 9._ 729 \ 
64 3 -4-^5 1 12 ' 12' 12' 12' 8' 2' 8' 8' 8' 2621442 B J 

45 v^sinCf) p /13 17 25 29.1 5 3 9 5 . 729 \ 
512 Tjfwl) 4 5 ^24' 24' 24' 24' 2' 8' 4' 8' 4' 2621442 s ' 

4_A 1 7? f' 2 . 5 1 Z 4. 5 3 7 9 5 11 . _ 729 \ 
~Tl6 F^ 5 ^ 6 ^' 6' x ' 6' 3> 8' 4' 8' 8' 4' 8' 2621442 s > 

693 r (4) sin (| ) p /19 23 31 35.3 7 5 11 3 . 729 \ 
8192 4^5^24' 24' 24' 24' 4' 8' 4' 8 ' 2' 2621442 6 > 

63 %/2 p /ll 13 17 19.7 9 11 3 13. 729 \ 
4096 8 r 4 5 U2' 12' 12' 12' 8' 8' 8 ' 2' 8 ' 262144^^ 

I 1989 v^sinCir) p/25 29 37 41.9 5 3 13 7._ 729 \ 
~ 1 ~262144 ^ r( 3, 4-^5 V 24' 24' 24' 24' 8' 4' 2' 8 ' 4' 2621442 s / 
^ 1 v 4 1 

4 16 r (g) sin (~r) p/I 1 » 29 II 1-1 1 2 1 1 I 9 . 16384 \ 

5 25 I 6^7^40; 20' 40' 40' 20' 40' 5' 10' 5' 2' 10' 5' 10' 97656252 s J 

48 r (5) sin (f) p/13 1233319 43,_3_2134_9 11 . 16384 \ 
125 I 6-^71.40' 20' 40' 40' 20' 40' 10' 5' 2' 5' 5' 10' 10' 97656252 s / 

^2 -> 7T 

112 csc(^) sin(| ) ^^ 11273721 47. 2131111 6 . 16384 
'625 3 



5 ; "'"i 5 ; p f!I 112T372147.213_/_ill6. 16384 \ 

12 .3 " 6-^7140' 20' 40' 40' 20' 40 ' 5 ' 2 ' 5 ' 10' 10' 10' 5' 97656252 s > 
z & n g ; 



I 176 csc (f ) sin (ir) p/211JJ3J,4123 52.13_7_4116 13. 16384 \ 

'3125 ¥^7T\ 6 ' r7 U0' 20' 40' 40' 20' 40 ' 2 ' 5 ' 10' 5' 10' 5' 10' 97656252 s ' 

z & 

6384 £(|)s™(xj 7? /29 17 39 49 27 59. J7_ 4 _9_ 6 13 7 3. 16384 \ 
78125 M 6^7^40) 20 ' 40 ' 40 ' 20 ' 40 ' 10 ' 5 ' 10 ' 5 ' 10 ' 5 ' 2 ' 97656252 s > 

14352 r (75) sin (f ) p/B19«5329 63.4_9.111373 8 . 16384 \ 
390625 W 6^7^40; 20 ' 40 ' 40 ' 20 ' 40 ' 5 ' 10 ' 10 ' 10 ' 5 ' 2 ' 5 ' 97656252 s / 

26928 c^c(-f : )sin(f ) p /37 21 47 57 31 67._9_ 11 6 7 3 8 17. 16384 
1953125 42 3 6^7^40) 20 ' 40 ' 40 ' 20 ' 40 ' 10 ' 10 ' 5 ' 5 ' 2 ' 5 ' 10 ' 97656252 s 



2 



35464 csc (f ) sin (-Tf) p /41 23 51 61 33 71.11 6 13 3 8 17 9. 16384 - 
9765625 W77i\ 6 7 U0' 20 ' 40 ' 40 ' 20 ' 40 ' 10 ' 5 ' 10 ' 2 ' 5 ' 10 ' 5 ' 97656252 s - 



»» T($) 
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Table 2: The Fox's //-function for 1 < a < 2. 



a H{z) = azf(z) 



4 6ljl csc(£)csc(^) r(|f)r(j)r(f) F , 7 ll 19 23. 1 3 1 3 7. 729* 8 \ 

3 ^ I 2 csc(^) r(g)r(i|)r(|) 4J ^24' 24' 24' 24' 4' 8' 2' 4' 8' 262144^ 

ji ^(^Wff ) r(||)r(|i)r(|) p/13172529.1539 5. 729 Z 8 \ 

384 Z csc(|) r(^)r(i^)r(^) 4 5 ^24' 24' 24' 24' 2' 8' 4' 8' 4' 262144 J 

I 77 „4 gf£(fi Mf ) rTi)^)^!) pf 19 23 31 35. 3 r 5 11 3. 729*8 n 
~ 1 ~1024 Z csc(^f) r(g)r(||)r(|) 4 5 ^24' 24' 24' 24' 4' 8' 4' 8 ' 2' 262144/ 

221^ „6 csc(^)csc(|f ) r(||)r(||)r(|) F /25 » 37 41. 9 5 3 13 7. 729z s x | 
196608 2 csc(|) r(||)r(||)r(|) 4 5 ^24' 24' 24' 24' 8' 4' 2' 8 ' 4' 262144 J I 

2 h 7 ^(|M|Wf) r(f )r(^)r(f )r(|) , 9 7 19 29 17 39. 1 3 21 7 4 9. 16384^% 

4 v^f^l CSC (|) CSC (1|) CSC (^L) r(^)r(§)r(|§)r(§) 6^40' 20' 40' 40' 20' 40' 5' 10> 5' 2' 10' 5' 10' 9765625 J 

72§ r 2 ™=(ffiW&W|*) r(§)r(f )r(f )r(f) p ,17 11 27 37 21 47.2 1 3 J_ _9_ 11 6. _ 16384*^ 
100 z CSC (^L) CSC (JL) CSC (|) r(f )r(§)r(§)r(^) 6 7 U0' 20 ' 40 ' 40 ' 20 > 40 > 5' 2' 5' 10 ' 10 ' 10 ' 5' 9765625 ^ 

_i 1 _^4 csc (f ) csc (ir) 77 (5 3 7 1 9 5 11.3 X 4 J. 11 6 13 7 . 16384z 10 \ 

2 $j csc (JL) csc (|) csc (^I) csc (^I)7 r 8l 8 , 4 , 8 ' - 1 ' 8' 4' 8 ' 5' 10' 5' 10' 10' 5' 10' 5' 9765625 J 

2§897 r 6 " c (^W^Wlg) Jr (f)r(i)r(f )r(^) p ^ 19 43 53 29 63.4 I U 13 7 3 8 . 16384z 1 ° v 
31250 z CSC (^) CSC (1|) CSC (^L) r(§)r(§)r(li)r(f) 6 r 7l 40 , 2 ' 40 ' 40 ' 20 ' 40 ' 5' 10 > 10 ' 10 ' 5' 2' 5' 9765625 > 

, 2*4433 „8 csc (^) csc (|g) csc (w : ) £(iHmHljHAj F f 4 ! 23 51 61 M 7i.ll 6 13 3 8 17 9. 16384z 10 - 
3125000 2 csc (JL) csc (|) csc (3l) r(§)r(^)r(§)r(|) 6r7 Uo' 2 0' 40 ' 40> 20 ' 40 > 10 > 5' 10 > 2> 5' 10 > 5' 9765625 



2MT7 1 Wlf)^(lf) r (ln) r (f ) r (j) r (M) r (A) 



IX 



11 17 23 29 41 47 53 59.1 1 1 _5_ 1 2 3 5 11 . 9765625-z 12 \ 
60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 6 ' 4 ' 3 ' 12 ' 2 ' 3 ' 4 ' 6 ' 12 ' 8707129344 > 



5§2lh ^2 ^(jfW^csc^jcsc^) PC1 11113171921 23.1 _5_ 1 7_ 2 5 11 13 7. 9765625z 12 \ 
320 ^ CSC (52L) CSC (JL) 8^9^20' 20' 20' 20' 20' 20' 20' 20' 3' 12' 2' 12' 3' 6' 12' 12' 6' 8707129344 ^ 

10*1729 r 4 ^(g5wb) C 4%W^)r(i)r(i)r(g)r(|)r(i|) 
+ 2488320 z r(^)r(f)r(S)r(|i)r(^^ x 

7? I'M 37 43 49 61 67 73 79.1 7_ 2 3 5 13 7 5 4. 9765625z 12 \ 
A 8 r 9V 60' 60' 60' 60' 60' 60' 60' 60' 2' 12' 3' 4' 6' 12' 6' 4' 3' 8707129344 J 

124729 j ^(f)^^)^(t)^(lf) ^§) r (f ) r (f ) r (f ) r (A) w 
8957952 csc(!f) r( f )r(|)r(|)r(|)r(M J x 

'41 47 53 59 71 77 83 89. 2 3 5 11 7 5 4 17 3 . jm^^^ 2 \ 



X 8-^9 160' 60' 60' 60' 60' 60' 60' 60' 3' 4' 6' 12' 6' 4' 3' 12' 2' 8707129344 ^ 

55 2*429 8 csc (fl ) csc (M) cgc (if ) csc (ff ) F/17 192J_2^272^313^.5U137417319 5 . 9765625z 12 ' 
655360 ■* CSC (JL) CSC (&2L) 8-^9^20) 20 ' 20 ' 20 ' 20 ' 20 ' 20 ' 20 ' 6 ' 12 ' 12 ' 6 ' 3 ' 12 ' 2 ' 12 ' 3 ' 8707129344- 

10*596932063 i e.c(l)c,c(|w^M^)r(ii)r(|)r(|)r(|)r(||) , 

16717688340480 Z csc{^) TMW&&&i^ ' 



61 67 73 79 91 97 103 109.13 7 5 4 3 19 5 7 11 . 9765625z 12 ' 
60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 60 ' 12 ' 6 ' 4 ' 3 ' 2 ' 12 ' 3 ' 4 ' 6 ' 8707129344 - 
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8 Conclusions 



In this article we presented analytic expressions for the a-stable p.d.f.'s for rational values of 
the index a G (0, 2]. These p.d.f.'s can be viewed as solutions of the spatial anomalous diffusion 
equation subjected to a Dirac delta initial condition. We established their connection to the 
Fox's if- function for the most general law S a (/3, K a , r). The characteristic function of the a = 1 
stable law was also reproduced by solving a suitably chosen fractional diffusion equation. An 
alternative way of classification which captures the infinite divisible character of stable laws 
was proposed. The rationality of the index allows us to write closed expressions for the p.d.f.'s 
in terms of generalized hypergeometric functions. This method recovers known results, such 
as the Cauchy p.d.f. for a = 1, the Holtsmark p.d.f. for a — | and the normal p.d.f. for 
a = 2. When a takes an arbitrary rational value in (0,2] new p.d.f.'s are derived generalizing 
and unifying previous results. 



Appendix A 

The proof of (31) is straightforward provided that we first show that 

(i) If / and g belong to the space AC^lQR) 12 with boundary conditions limuuoo f^ k \x) = 
= lim^i^oo g [h) (x), k = 0, ■ • • , [a] - 1 then 

UWCaO) 9(x)dx = [ f(x) ( x V^g(x)) dx (A.l) 
i J i 

where I = (— oo, oo). 
Proof of (A.l) 

\ ^7 r ^Srdu) g(x)dx 



j °° " L\[a\ + 1 - a) Jx \axw~ rJ - j u 



1 


r(M 


+ i 


— a) 


(- 


!)[«] 


+i 


r(M 


+ 1 


— a) 


(- 


.!)[«] 


+i 


r(H 


+ 1 


— a) 



fix - u) J \ d^ +1 

n—du , r r-: g(x)dx 

J dx^ +1 

t/ . / d^ +1 f°° g{x) , \ . 
m [d^J s (x-s)^ dX ) dS 

f(x) ( x V^g(x))dx. (A.2) 



l 



(ii) The two identities hold 

-ooV*e ipx = {ip) a e ipx = |p| Q e^ sign(p) e ipx , 
x V^e ipx = (-ip) a e ipx = \ p \» e -^sn(p) e ipx_ 



Proof of (A. 3) 
Using (30) we have 

V a J px = — ^ r -^rre ipx / — -^du 







a—\a\ 



r([a] + l -a) \dxW +1 
(ip) a e ipx , p>0. (A.4) 



12 This space consists of all functions / which have continuous derivatives up to order [a] - 1 on 1 with 
/(H-i)(>) 6 AC(R). We also recall that AC{R) C 5(M). 
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In the derivation of (A. 4) we have used the integral formulas [25] 

f°° m-1 / \j r (^) fl™\ 
I ar cos ra; ax = cos — 

Jo V v 2 ; 

/ x^ 1 sin(px)dx = sin ( — \ , for p > and < \Re/i\ < 1. (A.5) 

Jo p** V 2 / 



Appendix B 

We first prove the identities 
(1) 



lim ( \x\ e - -5(x)] = (B.l) 
->o+ V e / M 



Proof of (B.l) 

Differentiating with respect to x the expression 



\x\ e — r 

ln(|x|)sign(x) = lim ( ) sign(x) (B.2) 



and taking into account that dsi ^ x ^ = 25(x) we easily derive (B.l). 

(2) 

= (27Tbl)^ r(1 + ^ ^ + ^ ' a>_1 ' (B ' 3) 

Proof of (B.3) 

This is a direct consequence of (A.5). 
The Fourier transform of ln(|x|) is 

F[H\A)\ = lim^f j-i^ UhmfM, 2 r(1 _ e)! infl5) 



-55f + C7 7*(P) ( B - 4 ) 



where C 1 is the Euler-Mascheroni constant. In the derivation of (B.4) we made use of the 
r(l + z), \z\ < 1 [27] expansion 



oo 



r(l + z) = ^a n z n (B.5) 

n=0 

where the coefficients are given by 

n 

a = 1, na n = -7<2 n _i + ^(-l) fc a n -fcC(*0- (B.6) 

fc=2 
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